Tavis-Cummings model and collective multi-qubit entanglement in trapped ions 



O 

o 

(N 
£ 

(N 
(N ' 

> : 

00 . 

^t- ■ 
o ■ 
in : 
o . 
\o ■ 
o ■ 

^: 

9 L,: 
^— > ■ 

G ■ 

g : 
cr 



A. Retzker, 1,2 'B E. Solano t, 3 < 4 and B. Rcznik 5 

1 Institute for Mathematical Sciences, Imperial College London, SW7 2PE, UK 
2 Q0LS, The Blackett Laboratory, Imperial College London, Prince Consort Rd., SW7 2BW, UK 
3 Max-Planck-Institut fur Quantenoptik, Hans-Kopfermann-Strasse 1, D-85748 Garching, Germany 
'*Secci6n Fisica, Departamento de Ciencias, Pontificia Universidad Catolica del Peru, Apartado Postal 1761, Lima, Peru 
5 Department of Physics and Astronomy, Tel-Aviv University, Tel Aviv 69978, Israel 

We present a method of generating collective multi-qubit entanglement via global addressing of 
an ion chain performing blue and red Tavis-Cummings interactions, where several qubits are coupled 
to a collective motional mode. We show that a wide family of Dicke states and irradiant states can 
be generated by single global laser pulses, unitarily or helped with suitable postselection techniques. 
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I. INTRODUCTION 

Multi-partite entangled states play an important role 
in quantum information. They are useful in various quan- 
tum information applications, such as in Heinsenberg- 
limited spectroscopy [lj, secure communication [2]], and 
various schemes related to "one-way" quantum comput- 
ing Fresh theoretical developments on the genera- 
tion of multipartite entangled states show that sequen- 
tial techniques may prove to be general and practical for 
building arbitrary multi-qubit states For instance, 
a recent experiment [5] has realized a W state of eight 
qubits, encoded in the internal ionic levels, by perform- 
ing a sequence of two-qubit gates on different ion pairs. 
However, given a set of available interactions in a physical 
system, there are particular families of entangled states 
that could be built globally and in fewer steps [6| . In the 
context of cavity QED (CQED), for example, the cou- 
pling of a single cavity mode with a two-level atom, the 
Jaynes-Cummings model (JC), can be extended to the N 
atom case, leading to the Tavis-Cummings model, with 
different dynamics and entanglement features 0, Q • 

In this article, we study methods of generating specific 
classes of multi-qubit entangled states in trapped ions 
with collective interactions, which are potentially faster 
and more efficient than individual techniques. They con- 
sist of two key ingredients: firstly, the use of global rather 
than individual addressing of ions and, secondly, the pres- 
ence of invariant subspaces, i.e., combined (vibronic) in- 
ternal and motional finite subspaces that are closed with 
respect to certain dynamical operations. 

In Sec. [ill we describe realistic collective vibronic in- 



teractions coupling the internal degrees of freedom of ./V 
ions with a collective motional mode. Specifically, we 
consider the blue and red excitation versions of the Tavis- 
Cummings model, taking distance from usual predictions 
in the Dicke model. In Sec. IIII1 we study the invariant 
subspaces, associated with the proposed interactions, in 
the search of classes of multipartite entangled states that 
may be efficiently generated. It will turn out that one of 
them is the family of symmetric Dicke states [i| [lo|, II j , 
from which the W state is just a one-excitation particu- 
lar case. In Sec. IIV1 we consider the family of entangled 
states that could be generated by means of purely uni- 
tary global operations and, in Sec.[V] the ones that could 
be generated by using postselection. 



II. COLLECTIVE MAPS 

Let us consider N ions in a linear Paul trap, cooled 
down to their collective motional ground state. We will 
not concentrate on a specific experimental setup , and 
our derivations could be applied to any ion-trap device. 
The free-energy Hamiltonian, Hq, describing the N two- 
level ions and their motion around their equilibrium po- 
sitions is 
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Here, cr™ are z-components of Pauli spin vectors describ- 
ing the two-levels with energy gap luq, while aj and at 
are the annihilation and creation operators for the nor- 
mal modes with frequency Vj. The interaction between 
the internal degrees of freedom of each ion and a col- 
lective motional mode can be induced by laser light of 
frequency u>, yielding [13 1 



Hfa t — h\ n cr™ cos(kx n 



LOt + 4> n ). 



(2) 



2 



Here, A„ is the coupling strength between the laser and 
the n-th ion, cr™ are x-components of Pauli vectors, k is 
the laser wave vector, x n is the displacement operator 
with respect to the equilibrium position, and <f> n is the 
phase of the laser at the location of the n-th ion. 

We will study the case of homogeneous laser excitation, 
A n = A, Vn, and of near resonant coupling, u> ~ ujq. For 
the sake of simplicity, we will also consider all (j) n = 0, 
although this may play an important role when mak- 
ing experimental considerations. In this case, the Hamil- 
tonian in the interaction picture, after a rotating-wave- 
approximation RWA) with respect to the two internal 
levels, reads 13[ 14 [ 



H l = ^£«e-^exp 
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a j e- iu i t )) + R.c. 



(3) 



where b n j denote the amplitudes of the j-th normal mode 



of the ion chain in the position expansions, 8 = uj — loq, 
and m is the ion mass. In the Lamb-Dicke limit, where 
all Lamb-Dicke parameters rjj — 



exponential can be expanded and set for a RWA with re- 
spect to the phonon field. In that case, when the laser fre- 
quency is tuned to a particular collective motional side- 
band frequency, u) = ujq ± Uj, we obtain blue and red 
sideband transition Hamiltonians 

hL* = Jaj£M°+ } + *- a i)> ( 4 ) 

n 

H Ld = ^£M*+«* + *-«5)> (5) 
n 

where Xj = r/jX/2. The interaction of Eq. (O appears 
naturally in the context of CQED, where a bunch of 
atoms interact inhomogeneously with a cavity mode and 
the counter-rotating terms are neglected in the RWA. 
The dynamics in Eq. ((4]) is not usual in CQED but can 
be easily engineered in trapped ions. Only when Vj corre- 
sponds to the center-of-mass (COM) mode frequency v\, 
we have b n \ = b\ , and we can define the collective angular 
momentum operators L + = J2 n <T + an d L- — J2n a - ■ ^ n 
that case, we could rewrite the Hamiltonians of Eqs. Q 
and ([5]) as 



^blue = fiAi(L + oT + L-a) 
H^ cd — hXi(L + a + L^aJ), 



(6) 
(7) 



where Ai = &1A1. The dynamics associated with the 
Hamiltonian of Eq. J7]) is named after Tavis and Cum- 
mings 0] , who developed the first analytical solutions for 
this model. When we consider a motional mode different 
from the center-of-mass one, we could always define L\ = 



J2 n bnj&+ and L 3 _ = Yln^nj^-, but these operators do 
not satisfy the usual angular momentum algebra. If we 
define L z = a n z and L 2 (j) = L\ + \{L° + L J _ + lP_lP + ), 
with j = 0, 1, • • • N — 1, we get 



[L 

[L z ,L 2 (j)} = 
[L^L\k)]^Q 



(8) 



In fact, Lj_ can still be used to lower and raise the quan- 
tum numbers of L z , but they do not commute with 
L 2 {j). For the case of the center-of-mass mode, where 
all commutations relations are satisfied, we shall de- 
note the eigenstates of L 2 (l) and L z by \l,m), with 
I = N/2, N/2 - 1..., / > 0, and -I < m < I. States 
\l,m) are known as the Dicke states 



III. INVARIANT SUBSPACES 



Hamiltonian H^. cd conserves the total number of spin 
and phonon excitations, and commutes with the excita- 
are small, the tion number operator R = ^ m aj„ a m + L z + N/2, while 
Hamiltonian Hi-, conserves the difference between the 
spin and phonon excitations, hence, it commutes with 
B = Yl a m a m — L z + N/2. It is therefore possible to 
consider vibronic subspaces with a fixed number of ex- 
citations associated with R or B. If we concentrate on 
the case H = H' J rcd , we have the eigenstates |r, a) of R, 
where r = 0, 1,2- ••, and a denotes other degeneracy 
lifting quantum numbers. We then obtain the block di- 
agonal structure H J rcd = ©r=o°-^rcd( r )- ^ ne dynamical 
evolution that is generated by -ff;?ed leaves the subspaces 
invariant. 

We proceed to discuss certain examples of such invari- 
ant subspaces, for example, the one associated with the 
case j = 1. The smallest eigenvalue of R, r — 0, corre- 
sponds to the state H r =o = {\l = N/2,m= -N/2)\0)}, 
i.e., all atoms in their ground state and no phonons in 
the system. For the case r = 1, we have 



H r =i — Hi =N / 2 © 'Hl=N/2-l, 



(9) 



where 



H/=jv/2 = { 1^/2,-^/2)11), LV/2,-iV/2+l)|0) } (10) 
and 

n l=N/2 ^ = { \N/2 - 1, -N/2 +l,a = 1)|0), .., 

\N/2- l,-N/2+ l,a = N- 1)|0)}. (11) 

The quantum number a — 1, N — 1, lifts the (N — 1)- 
fold degeneracy of the states with / = —N/2 + 1. Hence, 
values of a enumerate the different angular momentum 
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multiplets. It is important to stress that H^ ed does 
not mix the different multiplets and, since L_\N/2 — 
I,— N/2 + l,a) = 0, there are no further transitions. 
This does not follow merely from the conservation of 
R, which does not forbid transition between the state 
\N/2 - 1, -N/2 + 1, a) |0), which has terms with one ex- 
cited atom, and a state with one excited phonon. This 
non-mixing property of the multiplets reflects the effect 
of quantum irradiance {{J EH , 16 1. The construction of 
higher r-number subspaces is straightforward. For in- 
stance, for r = 2 we shall have Ti. r =2 — T~(-i=n/2 ® 
7~Ii=n/2~i © 7~Li=n/2-2i etc. 

A key point in the present work is the use of subspaces 
which arc bidimensional. In this simple case, the evo- 
lution of the system resembles that of the well known 
Rabi oscillations. For example, let us consider the r = 1 



invariant subspace 7i r= 



l,l=N/2- 



We can start with the 



non-entangled state containing one phonon and with all 
the internal spins in their ground states. When we turn 
on the Hamiltonian H^ cd we obtain an oscillation between 
the states 



\N/2, -N/2) |1) » | N/2, -N/2 + 1) |0) 



(12) 



State \N/2, —N/2 + 1) is a symmetric combination of N 
terms, (| |U • • •) + | III •••) + ••• | • • • U1))/y/N, known 
as the W state. Similarly, we could make use of the in- 
variant space 7~t r =i,l=NI2—x,ot and, in that case, we would 
have the following oscillation 



| N/2, -N/2 + 1) |1) <-> \N/2, -N/2 + 2) 



(13) 



In the general case, the invariant subspaces can be of 
higher dimension, for instance if we start with n phonons 
in the multiplet I = N/2, the relevant states for r = n 
becomes, up to rotations induced by 7J ro d, 



-JV/2)|I)«->|-iV/2+l) 



-N/2- 



I) |0) . 
(14) 



So far, we have discussed invariant subspaces which 
are connected with the Dicke states and the collective 
angular momentum operators with j = 1. By tuning 
the laser to couple other motional collective modes, we 
can access other j subspaces. As we discuss in the next 
section, it is sometimes helpful to combine several steps, 
and in each step to couple a different phonon normal 
mode. For instance, we can start with the state that 
contains two different phonon excitations 



iu-i)ii)ji: 



(15) 



couple first the internal levels with the phonon in mode 
i and later with the phonon in mode j. This process 



connects us with the state L J + L l + \l — N/2,m — —N/2). 
It is useful to see that in this type of transitions we have 

L{_ L% \N/2,-N/2)=LlJ2 b m\ii-i T I - I) 

n 

= fcbmbnA |||... |) =n%|U ...|), (16) 

where in the last step we used the orthogonality of the 
normal modes. 



IV. DETERMINISTIC CREATION OF 
ENTANGLED STATES 

With the use of i? rc d, many relevant states can be cre- 
ated. We start with the state 



\N/2,-N/2)\l)^ =|||| ...1)11), 



(17) 



where the j-th mode is occupied by a single phonon, 
and the internal state is not entangled. A W-state 
\Wf) = \N/2,-N/2 + 1) can be created by applying 
a single collective 7r/2-pulse on the state of Eq. (fT7|) . 
This can be easily understood by recalling that the above 
initial state belongs to the bidimensional Hubert space, 
n r=1 , l=N/2 = {\N/2,-N/2)\l),\W?}\0}}. In principle, 
by precise control of the duration and intensity of the 
laser pulse, a W state can be created between a large 
number of ions. In fact, a IF-state shared by eight ions 
has been created recently using a multi-step sequential 
procedure based on individual ionic addressing [5J . In the 
present proposal, we would require the previous prepa- 
ration of a single phonon in the COM mode and the ap- 
plication of a single homogeneous global laser pulse. A 
related scheme in the context of quantum dots was dis- 



al. [17 1- 



cussed recently by Taylor et 

It is also possible to generate deterministically higher- 
excitation Dicke states using other bidimensional invari- 
ant subspaces. The r = 2 subspace 'H r =2.i=N/2-i.a is 
a two-dimensional space that is spanned by the states 
\N/2 - 1, -N/2 + 1)|1) - and \N/2 - 1, -N/2 + 
2) 1 0) ~ W2 . The first state above is equivalent, up to 
local operations, to the Wf = \N/2, -N/2 + 1), while the 
second state, contains terms with two excited atoms and 
is equivalent, up to local transformations, to the second 
Dicke state \N/2,-N/2 + 2) = . The construction 
of W2 can therefore proceed as follows. We first obtain 
as described above , using a single pulse. In the sec- 
ond step, we transform Wf -> \N/2 - 1, -N/2 + 1) by 
changing locally the phases of each ion. This step re- 
quires local addressing implementing local rotations. In 
the final step, we add a single phonon and apply again 
H^ ed to obtain W 2 N , the second member of the subspace 
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r = 2. Unfortunately, it seems that for higher-excitation 
Dicke states, e.g. , this "climbing the ladder" method 
requires also some interaction between the qubits. To 
overcome this difficulty we shall discuss other methods. 

We consider next extended examples of coupling to 
other modes and show that they can be used for generat- 
ing irradiant states 

[E EE EH- We start with the state 

tpo = \N/2, — iV/2)|l)„. involving one phonon in the j-th 
mode and all the internal levels in their ground state, 
then, we apply the Hamiltonian H 3 The conservation 
of R restricts the possible evolution to the subspace of 
states with r = 1, i.e., to the states L 3 + ij}Q and L^L^rpo, 
with k = 1,2,. ..TV. However, we notice from Eq. (fTB|) 
that only terms with k = j do not vanish, hence, the 
evolution leads to Rabi oscillations in the bidimensional 
Hilbert space {ipo, L J + ip }. In this way, we can generate 
the family of entangled irradiant states of the form 



Kill 



(18) 



Irradiant states are states that do not emit photons 
and are thus more robust to decoherence than radiant 
states. In our case, this property is due to the relation in 
Eq. (fl6|) . Since the coupling to the electromagnetic field 
is through the L± operators, as it is for the phonon field, 
the resultant states are irradiant For the the case 
of two spins, the resulting state is the EPR state. The 
experimental feasibility of irradiance and superradiance 
in ion traps was discussed and demonstrated by De Voe 
and Brewer 

Having produced certain irradiant states, we can use 
them as a starting point for the deterministic genera- 
tion of an additional class of states. Irradiant states in- 
troduce other bidimensional invariant subspaces. Since 
L_ \ipirr) — 0, the subspace {\ipirr} |1) ,L+ \ipi rr ) |0)} is 
an invariant subspace of Hamiltonian i? ro d and, therefore, 
the second state can be produced by Rabi flipping. This 
is a new kind of entangled state which is a superposition 
of states with two spins in the upper state, 



* oc ITI- IT II -I) + perm 



odd 



I IT IT II- l> + perm 



(19) 



For the case of four spins the outcome of this process 
is a GHZ state. First, we apply H^ cd , which couples 
the internal states with the higher collective mode, j = 
N — 1 = 3, and create the irradiant state 

l% III ... |> = ITIII) - IITII) + IUU) - HUT) • (20) 

In the next stage, we apply the j = red Hamiltonian 



and get 

L^iu...i> = iuu>-imT>, 

which is, up to a local operation, a GHZ state. 



(21) 



V. CREATION OF ENTANGLED STATES 
WITH POSTSELECTION 

In the previous section, we have discussed determinis- 
tic schemes for producing irradiant states as well as the 
lowest Dicke states (including the W state). However, the 
full family of Dicke states could not be generated using 
only collective unitary transformations. In the present 
section, we present another approach which is useful for 
producing the full set of Dicke states 



W k N = 



-1/2 



ITT ..Til- I) + perm 



\N/2,-N/2 + k) 



(22) 



The properties of the Dicke states may be of considerable 
interest in quantum information, and have been discussed 
recently by different authors [ill, HH, 0, 0- It 
can be shown that the von Neumann entanglement en- 
tropy, with respect to a bi-partite split of N qubits in 
a Dicke state, increases with k and saturates gradually 
for large k values. The behavior of the (mixed state) en- 
tanglement between two qubits 19( can be evaluated by 
considering, for example, the negativity which increases, 
almost linearly with k. 

The basic idea behind our approach is that while a col- 
lective unitary transformation cannot be used to create 
any Dicke state, a suitable choice of the initial phonon 
state can bring us very close to our goal. In this scheme, 
however, there will be always a small error due to mixing 
with other states. Therefore, unlike the previous exam- 
ples, we propose to postselect the phonon state in order 
to be certain that the desired Dicke state was produced. 

In order to create the Dicke state , we begin by 
preparing the initial state \N/2,—N/2)\k) Vo . We then 
apply the time evolution of the Hamiltonian H 3 which 
takes this state into the r = k invariant subspace. It 
turns out that at a certain time the probability distri- 
bution will be sharply peaked around a state with zero 
phonons and Wj? for the internal levels. By measuring 
the number of the phonons it is then possible to remove 
the admixture of with other states. A procedure to 
create and measure the number operator in an ion trap 
was introduced by different authors HE Elllii- Experi- 
mentally, motional Fock states were already produced in 
the lab [23j, although those techniques required a series 
of consecutive Rabi flips. 
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The entropy of the internal degrees of freedom 
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FIG. 1: Population of various states in Hilbert space as a 
function of time for 100 spins and 40 phonons. It can be seen 
that the population of the last state is maximal. 




FIG. 2: The entropy of the internal degrees of freedom as a 
function of time for 100 spins and 40 phonons. It can be seen 
that the final state is nearly pure (see Fig. 1). 
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FIG. 3: The probability of states 
a function of time. 
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The crucial ingredient in our proposed mechanism is 
that the purity of the state, prior to postselection, is 
high. The fact that the state containing zero phonons in 
Eq. (10) is produced with high probability is shown nu- 
merically below. The intuition behind this is that there 
is an analogy between these subspaces and the angular 
momentum subspaces of L 2 , though the commutation 
relations are different. The L x operator rotates the spin 
about the x axis, producing states with L z = ±1 with 
a probability one. In order to model this dynamics an 
analogy could be made between this dynamics in Hilbert 
space and the dynamics of a particle travelling between 
sites with different coupling strength. Since the couplings 
are higher at the middle and lower at the edges, the prob- 
abilities are maximal at the edges, see Figs. [TJ and [2j We 
therefore expect the Hamiltonian H rc( \ to rotate the state 
between the first and last state in Eq. (fT4)) with a prob- 
ability close to unity. 

We make now some further considerations concerning 
our scheme based on postselection. The only states which 
are created with high probability are \l z = m), where m 
is the number of phonons. This is due to the fact that 
the last state is created with high probability, therefore, 
the number of phonons in the first state determines the 
final state. The population of states in the subspace of 



Eq. (fT4|) starting with |— 2) |2) are shown in Fig. [3l where 
we observe that for specific times the desired state is ob- 
tained with high probability. A similar thing is observed 
in Fig. [4] starting with |— 5) |5). In spite of the fact that 
the number of excitations is not negligible compared to 
the number of spins, the purity of the final state is con- 
siderably high. This observation may prove very useful 
for generating Wj^ states. 

In order to increase the purity of the final state the 
number of phonons has to be measured and the vacuum 
state postselected. To achieve that goal, we consider a 
recently proposed technique for sorting a desired mo- 
tional Fock state \N) out of any motional distribution. 
This technique is based on a suitably designed vibronic 
scheme in a single ion, allowing for a restricted dynamics 
inside a chosen selected JC subspace { \g) \N+ 1) , | e) | N) }. 
To adapt it to our present work, we would need an addi- 
tional idle ion inside the chain, coupled to the motional 
mode of interest and specifically assigned to postselec- 
tion purposes. Together with the additional necessity of 
individual ion addressing for the sake of manipulation 
and measurement, these requirements for the idle ion 
are at reach by the state-of-the-art present technology in 
trapped ions [5j]. The proposed scheme described hith- 
erto can also be applied to create motional number states 
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in trapped ion systems. We have considered two schemes, 
one based on purely (deterministic) unitary operations 
and the other one based on an ulterior (probabilistic) 
postselection. Both schemes use the fact that the Tavis- 
Cummings model, in its blue- and red-excitation ver- 
sions, possesses invariant subspaces. In the deterministic 
case, the global laser pulses produce the desired entan- 
gled states after rotations in the associated bidimcnsional 
invariant subspaces. In the probabilistic case, the allowed 
rotations produce edge states that are very close to the 
desired entangled states, requiring a highly efficient post- 
selection technique. We believe that all proposed schemes 
are realistic and at within reach using present state-of- 
the-art technology in trapped ions. 
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FIG. 4: The population of the various terms. It can be seen 
that except for the first state, only the last state approaches a 
value close to I. The first state is a state with 10 spins down 
and 5 phonons 



via the Hamiltonian -ffbiuc, which will rotate the state in 
the proper subspace. Postselecting the spin state will 
yield the Fock state state \N) and the number of spins 
measured up would indicate the number of motional ex- 
citations N. 

VI. CONCLUSIONS 

In conclusion, we have presented methods of producing 
entangled states using homogenous global laser coupling 
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